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Abstract 

We present a method to calculate small stationary fluctuations around static solu- 
tions describing bound states in a (l + l)-dimensional Xcj)^^^ theory in a finite domain. 

D ' We also calculate explicitly fluctuations for the \(j) . These solutions are written in 

^ \ 

terms of Jacobi Elliptic functions and are obtained from both linear and nonlinear 

equations. For the linear case we get eingenvalues of a Lame type Equation and the 

> ■ 

Q^ ' nonlinear one relies on Hirota's Method. 

(N 

p. 

i> ■ 1 Introduction 

o 

oo . 

O ' It is well known that, for some applications, quantum fields can be thought as 

classical fields upon which are added quantum corrections [T] . In addition non-linear 

^ ■ field theories are nowadays powerfuU tools to describe fundamental physical theories 

c^ ' including cosmology, mainly for the infiationary model of the universe [2l [3] . On the 

other hand it is also well known that quantum systems when placed in finite domains 
(cavities) can present significant alterations on their behaviour. A step further is the 
study of interacting fields in finite domains [1]. A modern paradigm is the famous 
Casimir Effect [SI E] • In particular since these corrections are sensitive to boundary 
conditions an effect could expected on bound states of the physical system under 
consideration. 
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In a previous work Carrillo et al [7] obtained a number of solutions of a (1 + 1)- 
dimensional A(/)^-Theory in a finite domain in terms of Jacobi Elliptic functions. In 
this work we take a step further in the direction of a semi-clasical approach of the 
theory by studying small fluctuations around one of the above mentioned solutions 
describing its bound states. We think this is enough to exemplify our approach to this 
problem. Nevertheless in this work we do not pursue on quantization. Below we show 
these fluctuations are described by a non-linear equation. In section 2 we present the 
linear case and in section 3 the non-linear case which is solved by Hirota's method 
[8]. In section 4 we present some conclusions and perspectives for future work. 

We consider the Lagrangian density £ of a scalar field : 

c = ^d,d^<p + ^M'<p' - ^r^\ (1) 

where A is a coupling constant and fi = x,t. 

From ([1]) we derive the following equation of motion, 

- df,d''(j) + M^(t) - Xcj)"" = 0. (2) 

With the signature (H ), we can write the equation ((21) in dimension (1 + 1) as: 

^.^,M',-.r-0. ,3) 

The problem now is to solve the equation ([3]) for (p decomposed in a sum of a 
static solution (f>o{x) and a small fluctuation r]{x,t). In fact, in order to compare our 
results in linear case with the non-linear one, which will be presented in section 3 we 
can, without loss of generality, write 

0(x,t) = 0o(a;)+r/(x,t). (4) 

Using dl]) in ([3]) and taking the equations for (^o and rj we find: 



dx 



+ M'<Po - X<j)^ = 0, (5) 



and 



- I? + (Af ^ - ^-fn " - ^ E ('^j ro-'v' = 0. (6) 

Solutions for the equation ([5]) are, in general, given by Abelian Theta Functions 
using the First Integral formalism [9]. Now we want to solve the equation ([6]), and, 
for this we divide the analysis in two cases: 



2 The Linear Case 

Since r^ is a perturbation in the equation ([6]), we assume that r]{x,t) = eC,{x,t) 
where e is a small constant and then, the last term, which is of order O [rf) can be 
neglected. So the equation reduces to: 

Here we are interested in stationary solutions of form: 

i{x,t) = e'-'^{x). (8) 

So, Eq.(I7]) is reduced to equation for -0: 

^+(M^ + .^-An0nV^ = O. (9) 

Now, we define a parameter E by equation: 

^ = M^ + u\ (10) 

and then, the equation ([9]) takes the form: 



^'^ + f f - An^rM ^ = 0. (11) 



This equation can be interpreted as Schrodinger Equation for a potencial of the 
form 00 ~^. The problem now is, given a static solution 0o of the equation ([5]), find a 
general solution of the equation (ITT1) . 

On the other hand, A0^-Theory (ra = 3), the general static solutions of the classical 
equation of motion to the field (j) are given by sn-type elliptic functions [TU] 



, , MV2c (Mx I , \ , , 

00 X =± ^ sn ^^1 + ^/1^2^,/, 12 

where c is a parameter belonging to interval (0, |] and 

1 = ^ , (13) 

c 

From this relation, clearly / G (0, 1]. 

Now, substituting 0o(x) given by Eq. (fT2l) in Eq. fITT]) with n = 3 we get: 
^'^(«) /«, 2, n (1 + 



(ia^ 



6/sn^(a,/)- ^ ' E] i){a), (14) 



where we have used also Eq. (TTOj) and the change of variable 

Mx 



a 



In the literature the general form of Eq. (IT4l) is the Lame differential equation, 
which is given by [TT] 

iJ^ _ (m(m + 1) A;2 sn^fa, A;) + C) A(a) = 0, (15) 

where m is a positive real number, k"^ is the parameter of the Jacobian Elliptic Func- 
tion sn, and C is an arbitrary constant. Therefore, comparing (fT^ with (TT5|) we 
obtain, m = 2 and C = — -^r^ -E. 



The solutions of the Lame equation have been studied in the literature [12] . With a 
few algebraic manipulations we get the following eingenfunctions and their respectives 
eingenvalues: 



/i\ / \ /■ / 3 ,^ , , Mx ,, , Mx ,, 

(1) r]i{x,t) = eexp(2-y/^^Mt)sn(— ==,/)cn(-==,/) 

for ^l=(Y^^M^ (16) 



/ / 3/ , ^ , , Mx ,, , , Mx ,, 
(2) r/Kx,t) = eexp(z-^/^^Mt)sn(-^===,Z)dn(-^===,/) 

for col = (^) M^ (17) 

Mr Mr 

{3) vi{x,t) = ecn{^==j)dn{^=^,l) for ool = 0, (18) 

V 1 + t V 1 + t 



(4) r]i{x,t) = €exp{iw4t) I sn ( , /) 



2, Mx „ 1 + / + V/2 - / + 1 



VTTT ' 3/ 



for uj = { :^^ ) M^ (19) 



(5) 'r]i{x,t) = eexp(za;5t) sn (—==,/) - 



2, Mx ,, 1 + l-y/P -l + l 



VTTT 3/ 

2 / 1 + / + 2v//^"^TTT\ ,,2 

for 0.2 ^ ( _^_ j ^2_ (20) 



From the previous relations, it is worth to note that taking / = 1 in (jTSl) . ( 1191) 
and fITB]) . fITTI) . we recover respectively, the eigenfunctions and energy levels of the 
ground-state and the first excited state (static case) of the Dashen-Hasslacher-Neveu 
(DHN model) [13]. For the fifth case, we observe that for I = 1, uj"^ = 2M^, which 
corresponding to the continuum part of the espectrum. 



3 Bound States Fluctuations by Hirota's Method 

Now let us come back to the nonlinear model of the Equation ([6]). Again we will 
find the fluctuations 77, but now using a slight modification of the Hirota's Method 
[8]. This is done as follows. Firstly we make a dependent variable transformation, 
namely, 

rj(x,t) = ^4^. (21) 

In original Hirota's Method, a (3 parameter enters the above equation as a mul- 
tiplicative constant in order to simplify the differential equations obtained. In our 
method, this is not necessary, so the /3 parameter does not appear in the equation ( 12T|) . 

Thus, substituting (pTjl in ([6]), we get the following equation 
$'$xxx + (6 - 3n)<D<I>,<l>,, + (n - l)(n - 2)<l>3, - (n - l)(n - 2)<I>,<I>2+ 



{n - 2)$$,$ji + (2n - 4)$$^$,^ - <l>2$,„ + Af2$2$ 



- A*" V I 0^-'^$^$2^^-"'= = 0. (22) 

k=i Vv 

The Hirota's Method prescribes that $(x) can be expanded as a power series: 

<l>(x,t) = l+7/i+7V2 + ... , (23) 

where fi{x,t) are unknown functions {i = 1,2,3, ...) and 7 is an arbitrary constant 
[8]. Substituting fl23|l into fl22l) and equating the coefficients for every degree of 7 to 
zero, we obtain an infinite system of differential equations (we show below only the 
first two): 

/i,xxx - fi,.u + {M' - An0^i) A,, = 0, (24) 



and 



f2,xxx - f2,xtt + (M^ - Xtk/)"^-^) f2,x = 

= 2/l/l,,«-2(M2-An0^l)/l/l,,+ 
((3n - 6)/i,^^ - (n - 2)fi^tt) fi,x - (2ra - A)fi^th,xt - '2fi,xxxfi 

-A^^^r'/L. (25) 



For now, we especialize for A0^ theory. In this case the above equation reduces to: 
-$,«$' + 2$$^i$i + $$,$« - 2$,$^ + 

$^^^$2 - 3 $$,.$,.. + (2 - A)*^ + M2$,,$2 - 3A(^2$,<I>2 - 3A(/)o$^$ = 0. (26) 
The equations for /i and /2 are now: 

- fiMt + fi,cccccc + (M^ - 3 A (j)l) /i,, = 0, (27) 

+ (3/l,a:x — /l,tt) /l,a; " '^fl,tfl,xt — '^fl,xxxfl — '^^4'ofl,x (28) 

and so on for higher orders O (7^). 

Notice that the function fi will be determined by the previous functions {/i, /2, .../j_i} 
only. So, we have a infinite set of recurrent equations. The next step is the usual one, 
as the first part of this work linear case). Thus, defining the new function 

ax,t) = h,,, (29) 

the Eq. (1271) can be written as 

-6t + exx + (M2-3A0^)^^Q_ ^3Q) 

Observe that the above equation is pretty the same as linearized equation ([7]). 
Since ^ satisfies ([7]) we can then use the solutions (( TT6l) — ( l20l) ) to obtain the functions 
/i, which are given by: 



:i) /i(x,t) = -^exp(z^^^Mt)dn(-^=,0, 

for wl = (^\ M\ (31) 



(2) /i(x,t) = -^^exp(zy^Mt)cn(-^,/), 

for wl = (^) M^ (32) 



(3) /i(x) = ^^sn(-^,0 for wl = 0. (33) 



(4) Mx,t)=exp{t\ -^ Mt)! ^ 



i+l V IM 'VT+T 



for wl = i ^^ j M' . (34) 



(5) fi{x,t)=exp{i\ — — Mt) 



l + l \ IM ^^/TTl 

^, Mx ,,, {l + l-^/f-l + l)x 
Eiam-^==, D) 

l + l + 2VP-l + l\ ,.2 



for wi = (^ ^^ J M\ (35) 

where, in all cases the constants of integration are taken as equal to zero. 

The next step consists in the determination of the function /2. In this case, using 
the Eq. (1271) into Eq. (^ results the equation for /2: 

12,... - f2,.u + (m2 - SXcPl) h,. = {?>h,.. - h,u) k. - 2/i,t/i,.t - 3A(^2 j2^ (36) 

One can see that substituting the functions /i given by fl3T|) - fl35|) in above equa- 
tion, we obtain inonhomogeneous differential equations for the function /2. Unfortu- 
nately due to the its great complexity, it was not possible to obtain the contribution 
of /2 with this approach. Nevertheless we have obtained solutions for the nonlinear 
equation (jS]), considering only the contribution of the function /i. On the other hand, 
7 is an arbitrary constant in the Hirota's Method, since the series in equation (12^ 
is just a formal series [8]. On the other hand, as mentioned above, in Eq. (j4]), the 
field 77(2;, t) is a small fluctuation on a static solution (l)o{x). We have found solutions 
below whose amplitudes are of order 7 (flrst approximation). So, for our application, 
7 is a small number controlling the fluctuation amplitude. Using fl2T]) . we obtain the 
corresponding fluctuations 77, which are given by: 

7exp(z./^Mt) sn(-^, /) cn(-^, /) 

(1) m{x,t) = ^1 ^J^ ^1^±^, (37) 

l-7^exp(2J^Mt)dn(;^,/) 



7exp(z J^ Mt) sn(-^, /) dn(-^, /) 

(2) Vi{^,t) = y-^ ^^ ^^^±^, (3^ 

l-7^exp(7J^Mt)cn(^,/) 



7Cii(Ji|,,0dii(-M|.,() 



(4) ^K^,i) = / ^ ^^ ^ ^,(40) 



(5) r^,(x,t) = / ^ ^ '- ^m 

1 + 7exp(zy/i±^±^^f^Mt) ^£ezl±|pEEi) _ ^E(am^, /)J 

Now observe that, at first aproximation (0 < 7 << 1), the correspondent fluc- 
tuations of linear and non-linear cases coincides, which implies 7 = e in equations 
fl371— I 41j) . This show that linear case is, in fact, an first order approximation of the 
non- linear one. 



On the other hand, in [7j the authors showed that imposing Dirichlet boundary 
conditions on the field 0, given by (1121) . confined to a box of size L, it must satisfy 
the condition 



ML = AV1 + IK{1), (42) 

where 4:K{1) is a period of the Jacobi Elliptic Function sn(M, /) [15]. So, imposing the 
same boundary conditions for the fluctuations rj (same kind of confinement), implies 
that the same / = 1{L), obtained from above relation, must be used for r]. It is not 
difficult to see from our calculation, that, when r] is confined, the only consistent lin- 
ear and nonlinear fluctuations satisfying the condition (H2]) are those given by 016p . 
(fT7|l . fl37|) and fl38l) . for they have the Jacobi Elliptic function sn(w, I), from which the 
relation (142|) comes from. 

In the more general case, of the A^*^"*"^ theory, the first approximation for fluctu- 
ations can be calculated from Eq. (pT]) and Eq. (l2^ . It is given by 



r7(x,t)=7 ,, ,:^\„_2 - (43) 



:i + 7/1 

where the functions /i, obtained by Hirota's Method. 

4 Conclusions 

We developed an adaptation of the Hirota's Method in order to calculate small 
fluctuations around static bound states in a A(/)""'"^-theory. For the case n = 3, we 

8 



calculate those fluctuations in a Finite Domain (interval) in (l+l)-diniensions in 
order to solve a non-linear equation of motion. We showed that those fluctuations are 
written in terms of Jacobi Elliptic functions. Imposing that fluctuations must satisfy 
the same boundary conditions of the unperturbed bound states we can select a subset 
of solutions. We must stress that, in the limit / = 1 ones get the results by Knyazev 
on fluctuations around a kink solution [13], for the fluctuations in the Eqs. (1371) and 
(139|) . Even in this case (kink) we have obtained three new fluctuations, fl38|) . PUI) and 
(14T|) . Dirichlet conditions imply the only possibility for fluctuation (140!) is / = 1, so 
it does not exist in finite domains. The solution fHT]) there exists in any domain. 



Acknowledgments 

This work was partially supported by CAPES (Coordenagao de Aperfeigoamento 
de Pessoal de Nivel Superior). 

References 



[1] 

[2] 
[3] 
[4] 
[5] 
[6] 

[7 



[9 

[lo: 

[11 

[12: 

[13 



R. Jackiw, Rev. Mod. Phys. 49, 681 (1977). 

A. Linde - Particle Physics and Inflation - Kindle Edition - April 16 (2007); 

E.W. Kolb, M.S. Turner - The Early Universe- Addison Wesley - (1990); 

J.A. Espichan Carrillo an A. Maia Jr. Int. J. Theo. Phys. - 38(1999), 2185-2196; 

G. Plunien, B. Miiller and W. Greiner, Phys. Rep. 134, 87 (1986). 

V. Mostepanenko, N. N. Trunov, R. L. Znajek The Casimir Effect and Its Ap- 
plications - Oxford Science Publications Hardcover - April 10, (1997); 

J.A. Espichan Carrillo and A. Maia Jr. J. Phys. A: Math. Gen. 33 2081 (2000). 

M.J. Ablowitz and H. Segur, Solitons and the inverse scattering transform SIAM, 
Philadelphia, (1981). 

H.F. Baker, An Introduction to the Theory of Multiply Periodic Functions Lon- 
don, Cambridge University Press, (1907). 

J.A. Espichan Carrillo, A. Maia Jr. and V.M. Mostepanenko, Int. J. Mod. Phys. 
A (2000), 2645-2660; 

Z.X. Wang and D.R. Guo, Special Functions- Chap 11, World Scientific - 1989 

P.M. Arscott, Periodic Differential Equations Pergamon, Oxford, (1964). 

R. Dashen, B. Hasslacher and A. Neveu Phys. Rev. D 10, 4131, (1974). 



[14] M.A. Knyazev, Phys. Lett. A 153 110 (1991). 

[15] M. Abramowitz and LA. Stegun, Hanbook of Mathematical Functions Dover Pub- 
lications, INC., New York, (1972). 



10 



